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SUMMARY 


A method for determining directly the overall dynamic response in terms of 
bending moment of an aeroelastic model of a launch vehicle to random (buffet) 
aerodynamic forces has been applied to two model configurations in the Mach num- 
ber range from 0.50 to 1.20. The two configurations tested were a blunted cone- 
cylinder and a bulbous nose with a cylindrical afterbody. The cone-cylinder con- 
figuration was free from buffet over the range covered in this investigation. 
Significant buffet response was measured on the bulbous nose configuration. For 
this configuration the total buffet response was composed of components associated 
with the first and second free-free bending modes. 

Presented in the appendix is an analysis using the techniques of generalized 
harmonic analysis which develops in some detail the relationships necessary for 
predicting full-scale buffet bending moments from data obtained from wind-tunnel 
tests on dynamically scaled aeroelastic models. 


INTRODUCTION 


Buffeting of launch vehicles is an important aeroelastic problem area both 
in terms of the loads experienced by the entire vehicle and of the environment in 
which a vehicle component may be placed. In several recent launch vehicle fail- 
ures, buffeting has been suspected as a contributing factor. Since no adequate 
theoretical method exists for the determination of the magnitudes of the fluctu- 
ating loads on a particular configuration, the designer must depend heavily on 
experimental data. Wind-tunnel buffet studies on aircraft lifting surfaces pro- 
vided useful information in predicting the buffet characteristics on full-scale 
components. (See ref. 1, for example.) It is reasonable to assume that similar 
studies on launch-vehicle models might be equally beneficial. Consequently, an 
extensive wind-tunnel launch -vehicle buffet study program has been undertaken at 
the Amps and Langley Research Centers . Most of the data reported to date have 
been concerned with the fluctuating pressure input part of the problem. For 
example, in references 2 and 3 measurements have been made of the fluctuating 
pressures on several different bodies of revolution. In reference 4 fluctuating 
pressure data are presented for different size models of a large manned launch 


vehicle in two different test media. Reference k has provided some useful infor- 
mation on the scaling of fluctuating pressure data. Although pressure data play 
a significant role in understanding the overall buffet characteristics of a par- 
ticular configuration, application of these input data in the prediction of 
structural response is usually very difficult. Examples of such an application 
are given in references 5 and 6 . 

The present investigation was undertaken to develop a method for determining 
directly the dynamic response of an aeroelastic launch -vehicle model to random 
(buffet) aerodynamic forces and to develop relationships useful in predicting 
full-scale vehicle response from model test results. The method proposed is to 
conduct wind-tunnel tests on a dynamically scaled aeroelastic model of the launch 
vehicle. Such a model (acting as a mechanical analog) in a suitable wind tunnel 
which generates the correct aerodynamic input forces performs the very difficult 
time and space integrations producing the desired response which can be measured 
readily. 

This method has been applied to an aeroelastic model which was tested over 
the Mach number range from O.5O to 1 . 20 . The model response was measured in terms 
of bending moment by using a resistance-wire strain-gage bridge. The basic con- 
figuration tested was a slender blunted cone-cylinder designed to simulate some 
of the geometric and dynamic characteristics of a representative launch vehicle. 
The second configuration was a modification of the basic configuration in that 
the conical nose was replaced with a bulbous nose. The model was supported in 
such a manner that it was free to respond in simulated free-free bending modes. 

Presented in the appendix is an analysis, employing simple beam theory and 
the techniques of generalized harmonic analysis, which develops relationships 
useful in scaling model-response data to their corresponding full-scale values. 


SYMBOLS 


A(o>) square root of mechanical admittance 

a(t) nose amplitude 

C generalized damping coefficient 

C(o>), 0^(00), C] j (k) correlation functions 



aerodynamic damping coefficient 
critical value of damping 
control- system damping 


aerodynamic damping derivative, 
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structural damping coefficient 

time -dependent section lift coefficient 

time- dependent effective total lift coefficient 

natural frequency, 0)/ 2n 

mode shape based on unit nose deflection 

reduced frequency, Lx>/V 

total length 

effective moment arm 

generalized mass 

Mach number 

mass per unit length 

stagnation pressure 

dynamic pressure, ^pV 2 

reference radius 
time 

time interval 
velocity 

longitudinal coordinate 
total nose amplitude 
total weight 
angle of attack 
phase angle 

mean- square bending moment 
root-mean- square bending moment 
total root -me an -square bending moment 


total root-mean-square bending moment associated with response in first 
two free-free bending modes, + °2^ 

| nondimensional longitudinal coordinate, x/L 

| 0 any particular nondimensional longitudinal station 

®y(0,o)) power spectral density of nose deflection 

co circular frequency, 2jtf 

p air density 

p. mass ratio, MypjtR^L 

A product of lengths, 

Subscripts : 

F full scale 

M model 

m mth natural free-free bending mode, m = 1, 2, 3, . . . 

n nth natural free-free bending mode, n = 1, 2, 3, . . . 

Dots over symbols indicate derivative with respect to time. 

APPARATUS AND TECHNIQUE 

Wind Tunnel 

The Langley 8-foot transonic pressure tunnel was used in this investigation. 
This facility is a slotted-throat single-return wind tunnel capable of continuous 
operation throughout the transonic speed range at stagnation pressures from l/k 
to 2 atmospheres. Both test-section Mach number and density are continuously 
controllable . 

Model 

In order to examine the aeroelastic model approach in predicting full-scale 
buffet loads, a simple aeroelastic model was designed. The geometric and dynamic 
properties of the model were chosen to be representative of those scaled from a 
typical launch vehicle; however, for simplicity, the model was designed to be 
approximately a uniform beam as far as its mass and stiffness distributions are 
concerned . 


k 


Two different conf igurations 
were studied in this investigation 
and their geometries are shown in 
figure 1. Configuration 1 was a 
blunted cone- cylinder, having a 
14*5° semivertex angle conical 
nose mounted on a cylindrical 
afterbody with a fineness ratio 
of 8.99* Configuration 2 was a 
modification of configuration 1. 
The conical nose was replaced with 
a 29 ° semivertex angle conical 
bulbous nose (has a reflex angle 
downstream of the maximum diam- 
eter) . The afterbody fineness 
ratio for this configuration 
was 7.66. The ratio of maximum 
nose diameter to afterbody diam- 
eter was 1.6. 

Some of the details of model 
construction and mounting system 
are illustrated in figure 2. The 
shaker shown mounted on the sting 
in figure 2 was used for aero- 
dynamic damping measurements pre- 
sented in reference 7* The model 
"snubber" was used to restrain 
the model response when data were 
not being taken. To facilitate 
assembly, the model was made in 
three sections . The model con- 
sisted of a 0 . 03 - inch-thick fiber- 
glass cylindrical shell which was 
radially stabilized by several 
aluminum-alloy rings spaced at 
intervals of approximately 



(a) Configuration 1. 



(b) Configuration 2. 



(c) Longitudinal section through nose of 
configuration 2. 

Figure 1.- Line drawings of configurations 
tested. All dimensions are in inches 
unless otherwise noted. 


5 inches over the length of the 

model. The longitudinal bending stiffness of the model was governed by the stiff- 
ness of the fiber-glass shell. Scaled model natural bending frequencies approxi- 
mating those of a typical large launch vehicle were obtained by attaching lead 
ballast weights to the aluminum-alloy rings. The bulbous nose, which was made of 
wood and fiber glass, was slipped over the conical nose of configuration 1 to form 
configuration 2. (See fig. 1(b).) Some of the ballast weight was removed to com- 
pensate to some degree for the weight of the bulbous nose section. 


The model was sting mounted. It was attached to the sting by two pairs of 
soft flex springs. A photograph of a typical set of springs is shown in fig- 
ure 2(b). The springs were attached to the model near the nodal points of the 
first free- free bending mode in order to minimize the influence of the springs on 
the free-free modes. An effort was made to determine the effect of the springs 
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(a) General view of model configuration 1 and sting support. 


Attach 



(b) General view of a pair of model support springs. L-63-4738 
Figure 2.- Photographs of model and sting assembly, and support springs. 


to 


Attach to model .? 


on the frequencies of the first three modes of configuration 1 by supporting the 
model at the node lines on soft rubber supports. The natural frequencies and cor- 
responding node lines obtained from these tests when compared with those measured 
with the model mounted on the sting indicated that the effect of the attachment 
springs was negligible. 

The physical properties of the two configurations are presented in tables I 
and II and in figure 3* The model natural, frequencies are presented in table I 
along with the corresponding structural damping ratios. The structural damping 
was somewhat amplitude- and temperature-dependent. (See ref. 7*) The values 
presented in the table are average values. The mass properties are presented in 
table II. The mass distribution was obtained by weighing the model components 
prior to assembly. The generalized masses were calculated from the experimental 
mode shapes and mass distributions. Presented in figure 3 are the measured nor- 
malized mode shapes for the first two free-free bending modes for both configura- 
tions. No mode shapes are presented for the third free-free bending modes since 
these modes contained relatively large deflections in the "hoop 11 plane (involving 
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TABLE I.- MODEL FREQUENCY AND DAMPING CHARACTERISTICS 


Configuration 1 

Configuration 2 

Remarks 

f 

Cs/^cr 

f 

Cs/^cr 

11 

(a) 

12 

(a) 

Sting mode 

22 

(a) 

23 

(a) 

Rigid body 

pitching mode 

32 

(a) 

37 

(a) 1 

Rigid body 

translation mode 

85 

O.OO85 

79 

0.011 

First free-free 
bending mode 

223 

0.013 

206 

0.010 

Second free-free 
bending mode 

387 

(a) 

333 

(a) 

Third free-free 
bending mode 


a Not obtained. 


deformations of the cross section) . The mode shapes were determined by exciting 
the model at resonance at constant amplitude and measuring the response with a 
small velocity pickup. Also included in figure 3 are the first two natural free- 
free bending mode shapes and frequencies for a uniform beam having the same total 
mass and effective stiffness as the model configurations. As is shown on the 
figure, the mode shapes for the model tested in this investigation are very 
similar to those of a uniform free-free beam. 


Instrumentation and Calibration 

The instrumentation used in this investigation are shown schematically in 
figure L. The dynamic bending moments were indicated by a suitably calibrated 
four-active- arm resistance-wire strain-gage bridge bonded to the model fiber- 
glass shell. The location of the electrical center of the bridge is shown in 
figure 1 . The bridge output signal was amplified by a 3-kilocycle carrier ampli- 
fier. The amplified signal was recorded on a 6. 75 -kilo cycle frequency-modulated 
tape recorder and monitored on a true-root-mean-square vacuum tube voltmeter. 

An alternating- current calibration signal was also recorded periodically by using 
a variable-frequency audio oscillator. The relationship between strain-gage 
bridge output signal and applied bending moment was obtained from a static cali- 
bration. To verify that the strain-gage bridge was sensitive to bending strains 
only, two calibration procedures were used. With the model mounted as a simply 
supported beam, known shear loads were applied and produced known bending moments 
about the strain-gage bridge electrical center due to support reaction forces. 

The bridge output, or unbalance, was read by using a self-balancing potentiometer. 
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TABLE II.- MASS DISTRIBUTION 


Configuration 1 : W = 19.68 lb; 

M]_ = 0.108l slug; M2 = 0 . 074-0 slug 


Configuration 2: W = 21.29 lb; 

Mj^ = 0.1150 slug; M 2 = O.O562 slug 

1 

in, slugs/ft 

■ 


m, slugs/ft 

0 

.0136 

0.0155 


0 

.0180 

0.0180 

.0271 

.0155 


.0313 

• 04 l 4 

.O55O 

.0155 


.0446 

.0598 

.0818 

.0155 


.0720 

.0348 

.0958 

.4385 


.0983 

.0428 

•1355 

.0155 


.1120 

.4700 

.1542 

.1650 


.1511 

.0504 

.1682 

.8975 


.1695 

.2000 

.1720 

.0155 


.1832 

1.0362 

.2172 

.0155 


.1896 

.1424 

.2312 

.7485 


.2313 

.0453 

.2579 

.0155 


.2451 

.7761 

.2748 

.0155 


.2713 

.0416 

.3302 

.0155 


.2880 

.0296 

• 34-42 

.8375 


.3^23 

.0155 

.4121 

.0155 


.3560 

• 8375 

.4291 

.8695 


.4230 

.0155 

.5000 

.0155 


.4396 

.8695 

• 5140 

.8735 


.5090 

.0155 

.5845 

.0155 


.5230 

•8735 

.5985 

.8725 


.5925 

.0155 

.6700 

.0155 


.6060 

.8725 

.6870 

.8125 


.6765 

.0155 

.7475 

.0155 


.6930 

.8125 

.7620 

.1247 


.7520 

.0155 

.7800 

.6260 


.7665 

.1247 

.8460 

.0155 


.7840 

.6260 

.8550 

.0607 


.8490 

.0155 

.8870 

.0155 


.8580 

.0607 

.9000 

1.1275 


.8890 

.0155 

• 9460 

.0155 


.9020 

1.1275 

.9600 

.2257 


.9470 

.0155 

1.0000 

.0155 


.9610 

.2257 




1.0000 

.0155 
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A similar calibration was made where 
known couples were applied about the 
strain-gage bridge location. The 
results of this calibration were in 
good agreement with the one where the 
moments were produced by shear loads. 
Therefore, it was concluded that the 
strain-gage bridge was sensitive to 
bending strains only. 


Test Conditions and Procedure 


Itynamic bending-moment measure- 
ments were made on configuration 1 at 
angles of attack of 0° and k° through- 
out the Mach number range from 0.50 
to 1.20 at a stagnation pressure of 
1,800 pounds per square foot. Similar 
measurements at an angle of attack of 
0° were made on configuration 2 at a 
stagnation pressure of 600 pounds per 
square foot over the Mach number range 
from O. 5 O to I.I 5 . One data point for 
configuration 2 was taken at a Mach number of 0.90 and a stagnation pressure of 
1,800 pounds per square foot. Configuration 2 was not tested throughout the Mach 
number range at a stagnation pressure of 1,800 pounds per square foot since at 
this pressure the model response was so large that failure of the support springs 
resulted after only 1 or 2 minutes of exposure (snubbers retracted). The tunnel 
stagnation temperature was held constant at 120° F for all test conditions. The 
variation of the test dynamic pressure and Reynolds number for both configurations 
is shown in figure 5* The Reynolds number was based on the model diameter at the 
base of the conical nose section for both configurations. These lengths were 
5* 06 and 8.10 inches for configurations 1 and 2, respectively. 

The procedure for a typical test point was as follows: The test-section 

Mach number and density were adjusted until the desired flow conditions were 
obtained. The model "snubber" was then retracted, and the model was allowed to 
respond to the aerodynamic forces present. A ^5-second sample of the strain-gage 
output signal was recorded on the tape recorder. The root-mean- square value of 
the strain-gage bridge signal was read by using the true-root-mean-square vacuum- 
tube voltmeter. The model "snubber” was actuated and this procedure repeated 
until sufficient data had been determined to cover the Mach number range of 
interest . 


Data Reduction 

Some of the data recorded on the magnetic tape were reduced to power spectral 
densities by the use of an electronic analog analyzer. The tape-recorded data 
were analyzed in the frequency range from 0 to 5 00 cycles per second by using a 



Figure 4.- Schematic diagram of data- 
recording instrumentation. 
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6 # 79- c y c l e “P er ~ secon d bandpass filter. 
Overall root -me an- square values of the 
tape-recorded signals were obtained by 
an electronic analysis. Root-mean- 
square values obtained by using the 
electronic analysis were in good agree- 
ment with the root-mean- square values 
obtained by using the true- root-mean- 
square vacuum-tube voltmeter while the 
tests were being conducted. The elec- 
trical signals were converted to 
bending moment by using the strain-gage 
calibration constant. 

For configuration 2 the bending 
moment associated with response in each 
of the first two free-free bending 
modes was obtained by integrating the 
power spectra in the neighborhood of 
the resonant frequency of the desired 
mode. For the first mode the range of 
integration was from about 50 cycles 
per second to 125 cycles per second and 
was from about 160 cycles per second to 
2l+0 cycles per second for the second 
mode. Modal bending moments were not 
determined for configuration 1 because 
of the relatively low level response 
obtained for this model. 



Mach nusfcer 


Figure Variation of test dynamic pres- 
sure and Reynolds number with Mach 
number . 


ANALYSIS 


For wind-tunnel buffet studies on launch -vehicle models to be useful in 
predicting full-scale buffet loads, the necessary model— full-scale scaling rela- 
tionships must be known. A dynamic analysis of launch-vehicle buffeting has been 
considered in some detail in the appendix. In the analysis, based on the tech- 
niques of generalized harmonic analysis, the vehicle was assumed to be flying at 
constant altitude with a constant velocity. The only aerodynamic forces present 
in addition to the random component were damping forces proportional to the veloc- 
ity of the bending vibrations of the system. No loss of generality results from 
neglecting the aerodynamic inertia and spring forces since such forces usually 
are small when compared with their structural counterparts for a slender launch 
vehicle. (See, for instance, refs. 7 and 8.) Structurally, the vehicle was con- 
sidered to be a linear mult i-degree-of -freedom system. 

The final result obtained from this analysis for the total root-mean- square 
bending moment at some longitudinal station £ Q is 
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( 1 ) 


a T 2 (£ 0 ) = X 

n=l 


7 



n 


q 2 R 2 L 2 k C T 
V 



The total mean- square bending moment is a superposition of single-degree-of- 
freedom results, coupling terms having been neglected in the development and each 
mode being independently treated as a separate system. (See appendix.) The 
general term of the series in equation (l) is essentially the same solution that 
would be obtained for a single-degree-of-freedom system subjected to a random 
loading. (See ref. 1.) The right-hand side of equation (l) may be conveniently 

separated into three parts. The first part, the term l n 2 is the square of an 

effective moment arm. The second term, enclosed in brackets, is an a dmi ttance - 
type term. In particular, it is jt/2 times the maximum value of the mechanical 
admittance in the nth mode multiplied by the width of the admittance curve at the 
one -half -power point. The damping ratio which appears in the second term of the 
right-hand side of equation (l) is 


where 




+ 




aerodynamic damping ratio 



control- system damping ratio 



structural damping ratio 


In the notation of reference 7 the aerodynamic damping ratio is related to the 

effective aerodynamic damping derivative C* by 

h 



C cr 2pk 


The remaining terms are associated with the random aerodynamic loading. The func- 
tion C L n (k n ) is the correlation function of the random section lift coeffi- 
cients for the nth mode. Although not mathematically exact, a convenient way of 
thinking of this function is that it is the power spectrum of an effective random 
aerodynamic coefficient in the nth mode. 
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The use of equation (l) for scaling buffet loads is readily apparent. Since 

the reduced frequency k n = and C]^ n ( k n) would be the same for both a 

dynamically scaled aeroelastic model and the full-scale vehicle, the full-scale 
bending moment for the nth mode is related to the corresponding model value by 


CT 


2 



n, F 






( 2 ) 


The total mean- square bending-moment relationship between full-scale vehicle and 
model is 



Although the full-scale bending moment at a particular location along the 
vehicle can be determined directly from equation ( 3 ), the missile or launch- 
yehpcle designer needs to know the distribution of bending moments along the 
structure . Therefore, in making dynamic bending-moment measurements on a model 
using a single strain-gage bridge, a strain-gage location sensitivity factor must 
be determined since a bridge located say at the point of maximum bending moment 
in the first mode may not be very sensitive to moments produced by response in 
the second mode. The necessary sensitivity factors were determined for the first 
two modes of both model configurations by calculating the bending- moment distribu- 
tion due to inertia loading for motion in each mode. These distributions were 
calculated by using the experimentally determined mode shapes and mass distribu- 
tions and are presented in figure 6 for configuration 2 in terms of percent maxi- 
mum bending moment. The corresponding distributions for configuration 1 were very 
similar to those of configuration 2. The strain-gage bridge location sensitivity 
factors measured in terms of the ratio of bending moment about the strain-gage 
location to maximum moment in a particular mode are presented in table III. As 
is indicated by the data in the table, the strain-gage sensitivities were approx- 
imately the same for both model configurations. 


RESULTS AND DISCUSSION 


A method for determining directly the dynamic response of an aeroelastic 
launch-vehicle model to random (buffet) aerodynamic forces has been applied to 
two model configurations. The basic concept of this method is the direct meas- 
urement of bending moments induced on an aeroelastic model supported in such a 
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TABLE III.- STRAIN-GAGE LOCATION 



SENSITIVITY FACTORS 


(a) First free-free mode. 



Configuration 

Free-free 
bending mode 

St rain- gage location 
sensitivity factor 

1 

1 

0.436 

2 

0.794 

2 

1 

0.521 

2 

0.752 


manner that it is free to respond in its 
free-free bending modes. The results of this 
investigation are presented in figures 7 
to 11. 

The bending moments measured on the two 
configurations are compared in figure 7 a s a 
function of Mach number. Since tests were 
conducted at different levels of dynamic pres- 
sure, it was necessary to determine a param- 
eter which would remove the effects of dif- 
ferent flow conditions from the model data. 

A dynamic analysis of buffeting based on the 
techniques of generalized harmonic analysis 
(see appendix) indicates that for a system 
with only structural damping, the root-mean- 
square bending moment is directly proportional 
to dynamic pressure for a given Mach number. 
The results of reference 7 indicate that the 
aerodynamic damping for the two configurations 
studied in this investigation is small when 
compared with the structural damping for the 
range of flow conditions covered. Therefore, the aerodynamic damping was neglected 
and the bending-moment data are presented in the form of the ratio of root -mean- 
square bending moment to dynamic pressure, (it is of interest to note that in air- 
craft wing and/or tail buffet studies it is the structural damping that is usually 
assumed to be small. (See ref. 1.) For this case when the structural damping is 
negligible, the dynamic analysis of buffeting indicates that the root -me an- square 
buffet bending moments would vary linearly with the square root of dynamic pres- 
sure.) Data taken at angles of attack of 0° and k° are presented in figure 7 for 
configuration 1. The data for configuration 2 were taken at an angle of attack 
of 0°. As is indicated by the data in figure 7 there is quite a difference in 
level of bending moment when comparing the results for the two configurations. 

The bending moments obtained for configuration 1 are believed to be produced by 
the response of the model to residual wind-tunnel turbulence. This belief is 
substantiated by the agreement of the data for the two different angles of attack, 


(b) Second free- free mode. 

Figure 6.- Calculated bending- 
moment distribution for the 
first and second free-free 
bending modes of configuration 2. 





since it would be expected that buffet 
bending moments would vary with angle of 
attack. The response of the bulbous nose 
model is attributed to buffeting. It should 
be pointed out that the aerodynamic damping 
for configuration 1, although small when com- 
pared with the structural damping, was posi- 
tive (stable) over the test range covered. 

The aerodynamic damping for configuration 2, 
also small, does become negative (unstable) 
over a small Mach number range near 0.95* 

(See ref. 7.) The total damping, including 
both aerodynamic and structural components, 
was positive. An estimate of the total 
bending moment for configuration 2 which is 
actually due to buffeting can be made by 
using the data for the cone -cylinder model 
as a tare value since both models have simi- 
lar dynamic characteristics and approximately 
equal strain-gage location sensitivities. 

That is, that part of the total response for 
configuration 2 which is due to wind-tunnel 
turbulence would be approximately equal to 
the total bending -moment response measured 
for configuration 1. 

As pointed out in the section entitled 
"Analysis," it is necessary to determine the 
distribution of energy throughout the fre- 


Ijii r O Configuration 1, a = 0° 


□ Configuration 1, a = U° 



.5 .6 .7 .8 .9 1.0 1.1 1.2 


Mach number 

Figure 7*- Variation of ratio of 
total root-mean- square bending 
moment at strain- gage location to 
dynamic pressure vith Mach number 
for conf igurations 1 and 2 . 


quency spectrum in addition to determining 

the overall energy level if the model data are to be scaled to full-scale values. 
Presented in figure 8 are sample bending- moment power spectral densities for both 
of the model configurations investigated. These spectra are from data taken at 
a Mach number of 0.90, and a dynamic pressure of 603 pounds per square foot. As 
is seen from figure 8 there are peaks in the spectra associated with the natural 
frequencies of the two configurations. The maximum values of the peaks obtained 
for configuration 2 are considerably higher than those obtained for configura- 
tion 1. Most of the power is contained in the peaks associated with the free- 
free bending modes. In the frequency range from 0 to about ^0 cycles per second, 
there is a small amount of power. There are three resonant frequencies in this 
range, a sting bending mode, model rigid-body pitching mode, and model rigid-body 
translation mode. The sting and rigid-body translation modes are extraneous in 
that they would not be present for an actual launch vehicle in flight. A nonzero 
rigid-body pitch frequency would appear since aerodynamic and engine stiffnesses 
are present. Similar extraneous modes would appear for any wind-tunnel model. 
Since these extraneous modes do not have resonances near those of the free-free 
bending modes, and soft springs were used to support the model on the sting, the 
support spring effectively isolated the model response from sting motions. This 
result is illustrated in figure 9 which shows the variation with Mach number of 
the ratio of root-mean-square bending moment attributable to response in the 
first two free-free bending modes for configuration 2 to total root-mean- square 
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Frequency, cps 


Figure 8.- Sample power spectra. M = 0.90 and q = 603 pounds per square foot 





bending moment. The deviation of this 
ratio from unity is a measure of the 
energy contained in all modes excluding 
the first two bending modes. As is 
seen from figure these deviations 
are relatively small for all Mach num- 
bers except for M = 0.^0 • Since exam- 
ination of the power spectral densities 
indicated that most of the response 
other than that produced by the first 
and second bending modes was in the 0- 
bo 40- cycle-per- second range^ the dif- 
ference of the data from unity is due 
primarily to motion in the sting and 
two rigid-body modes. The data in fig- 
ure 9 indicate that the support system 
(including sting and flex springs) used 
for the model of this investigation was 
satisfactory in that no appreciable 
extraneous bending moments were 
produced. 
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A 2 + °2 
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Figure 9.- Comparison of root -me an- square 
bending moment at strain-gage location 
attributable to first two free-free 
bending modes with total root-mean- 
square bending moment at strain- gage 
location for configuration 2. 


For purposes of illusbra/tion fhe 
bending moments measured on configura- 

tion 2 at the strain-gage location have been extrapolated to full-scale values 
by using the scaling relationships presented in the Analysis section. The 
following assumptions have been made: first, a hypothetical full-scale vehic e 

has been selected such that the model -to-full-scale relationships are as shown 
in table IV: second, a typical orbital launch trajectory for such a vehicle has 
been assumed during which the vehicle weight, free-free bending frequencies, 
control- system damping, and aerodynamic damping vary with Mach number as indica e 
in table IV; finally, the full-scale structural damping was assumed to be the same 
as that measured for the model. The aerodynamic damping ratios for the first 
mode at Mach numbers from 0.90 to 1.2 were estimated from the aerodynamic damping 
measurements made on this configuration previously and reported m reference 7- 
Since no damping measurements were made for the first mode below - *9 ,.an 

none were made for the second mode at any Mach number, the aerodynamic damping 
below M = 0.90 in the first mode and throughout the Mach number range for t e 
second mode was assumed to be the same as that measured for configuration 1. 

( See ref . 7 • ) 


The predicted variation with Mach number of the full-scale root-mean- square 
bending moments at the strain-gage location for the first two bending modes is 
shown in figure 10. Also included in the figure is the variation of the total 
root-mean- square bending moment attributable to the first two bending modes. It 
may he seen from the figure that for a model test stagnation pressure of 
600 pounds per square foot in the Mach number range from about 0.»5 to 1.00 the 
second mode provides the major contribution to the total bending moment whereas 
outside this range the first mode predominates. Also shown m figure 10 are the 
predicted full-scale moments from the model test at M = 0.90 and 
p = 1,800 pounds per square foot. Although the amount of bending moment 
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TABLE IV.- ASSUMED FULL-SCALE CHARACTERISTICS 
AND SCALING RELATIONSHIPS 


Scaling relations assumed for M = 0.90: 


= 0.04 





= 6.714 X 10“ 5 


Vm 

Vp 


1.00 


Pm . 
— = 1.049 

pf 



25.0 


% 

ip 




1.049 


M 

lb/ft 2 

W, 

lb 

f l, 

cps 

f 2> 

cps 

fcl 

fcl 

fcl 

fcl 

fcl 

fcl 

0.50 

233 

361, 000 

2.78 

7.26 

0.0110 

0.0100 

0.0100 

0.0020 

0.0022 

0.0022 

.60 

310 

347, 000 

2.89 

7.54 

.0110 

.0100 

.0104 

.0024 

.0022 

.0022 

.70 

400 

336, 000 

2.98 

7.78 

.0110 

.0100 

.0108 

.0028 

.0023 

.0023 

.75 

450 

331,000 

3.03 

7.90 

.0110 

.0100 

.0110 

.0030 

.0024 

.0024 

.80 

500 

326, 000 

3.07 

8.01 

.0110 

.0100 

.0112 

.0032 

.0024 

.0024 

.85 

540 

322,000 

3.11 

8.11 

.0110 

.0100 

.0114 

.0034 

.0025 

.0025 

.90 

575 

317, 000 

3.16 

8.24 

.0110 

.0100 

.0116 

.0036 

.0070 

.0021 

•95 

6 08 

313, 000 

3.20 

8.35 

.0110 

.0100 

.0118 

.0038 

-.0008 

.0021 

1.00 

635 

310,000 

3.23 

8.42 

.0110 

.0100 

.0120 

.0040 

-.0002 

.0019 

1.15 

698 

299, 000 

3.35 

8.74 

.0110 

.0100 

.0126 

.0046 

.0005 

.0018 
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contributed by each of the inodes has 
changed slightly, the total bending 
moment is about the same as that deter- 
mined from the model tests at the lower 
stagnation pressure. It is seen that 
the peak value of the bending moment 
occurs at about M = 0.95* 

Of course, the designer is inter- 
ested in more than the variation of the 
buffet bending moment at some partic- 
ular vehicle station with Mach number. 
The distribution of the buffet bending 
moments along the vehicle must be known 
in order that the design bending -moment 
distribution (the buffet portion of 
which is usually based on a 3 a or 
other appropriate factor) may be deter- 
mined. In figure 11 the longitudinal 
distribution of the total full-scale 
root-mean- square bending moments due 
to response in the first and second 
free-free bending modes are presented 
for Mach numbers of 0.70 and 0-95 • 

These data were generated by using the 
results for the individual modes pre- 
sented in figure 10 in conjunction with 
the bending-moment distributions given 
in figure 6. As is seen from fig- 
ure 11, the two distributions are some- 
what dissimilar. At the lower Mach 



Mach number 


Figure 10.- Variation of full-scale root- 
mean- square bending moment at strain- 
gage location with Mach number obtained 
from tests on model configuration 2. 


number the major contribution to the 
total moment was made by the first 

mode. As a result, the shape of the moment distribution was very similar to the 
first-mode moment distribution presented in figure 6. At M = 0.95 the contri- 
bution of the second mode was substantial, and this condition resulted in the 
"saddleback" shape of the moment distribution. This figure illustrates the 
importance of higher bending modes in buffet bending-moment measurements on 


launch vehicles . 


CONCLUDING REMARKS 


A method for determining directly the overall dynamic response in terms of 
bending moment of an aeroelastic model of a launch vehicle to random (buffet) 
aerodynamic forces has been applied to two model configurations in the Mach num- 
ber range from 0-50 to 1.20. The two configurations tested were a blunted cone- 
cylinder and a bulbous nose with a cylindrical afterbody. The cone-cylinder 
configuration was free from buffet over the range covered in this investigation. 
Significant buffet response was measured on the bulbous nose configuration. For 
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Figure 11.- Distribution of full-scale root-mean- square bending moment attributable to first 
two free-free bending modes at Mach numbers 0.70 and 0-95 obtained from model tests on 
configuration 2. 




this configuration the total buffet response was composed of components asso 
ciated with the first and second free-free bending modes. 


Presented in the appendix is an analysis using the techniques of generalized 
harmonic analysis which develops in some detail the relationships necessary for 
predicting full-scale buffet bending moments from data obtained from wind-tunne 
tests on dynamically scaled aeroelastic models. 


It remains desirable, of course, to evaluate the aeroelastic model 
by direct comparison of* model and full-scale result sj however, sui a e 
data for this purpose are not available. 


approach 

flight 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 

Va., August 8 , 1963. 
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APPENDIX 


DYNAMIC ANALYSIS OF BUFFETING WITH EMPHASIS ON 
SCALING MODEL DATA TO FULL-SCALE RESPONSE 


The purpose of this appendix is to examine a launch vehicle which is sub- 
jected to a random (buffet) aerodynamic loading. The method of generalized har- 
monic analysis, which is treated in some length in reference 9 and was first 
applied to the analysis of buffeting in reference 10, provides the basis for this 
study. The vehicle is assumed to be flying at constant altitude with constant 
velocity. The only aerodynamic forces present in addition to the random component 
are damping forces proportional to the velocity of the bending vibrations of the 
system. No loss of generality results from neglecting the aerodynamic inertia 
and spring forces since such forces usually are small when compared with their 
structural counterparts for a slender launch vehicle. 

The set of differential equations which govern the free-free bending vibra- 
tion characteristics of the system under consideration is 

M n *d n (t) + Cnajt) + a»n 2 M n a n (t) = qRL f c L (|,t) h n (|) d£ (n = 1, 2, 3, . . .) 

J 0 

( 1 ) 


where 

&n(t) nose deflection in nth bending mode 

C n generalized damping coefficient in nth bending mode, including aero- 

dynamic, control system, and structural components 

M n generalized mass in nth bending mode 

CDn natural circular frequency in nth bending mode 

The right-hand side of equation (l) is the generalized random aerodynamic load 
expressed in coefficient form. The function c L (£,t) is the random section 
lift coefficient and R, L, q, and h n (|) are, respectively, reference radius, 
vehicle length, free-stream dynamic pressure, and mode shape of nth bending mode 
referred to unity at the vehicle nose. The solution of equation (l) for the 
nose deflection is 
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a- 


•n(t) 


qRL / 

J 0 


c Ti (|,t) h n (|) d| 


»w 


i - (2- T 

o>n ) 


+ i 


2 0(c c r) n 


( 2 ) 


or 


a n (t) = 


qRIA n (a)) c L n (t) e 


-10 


n 




(3) 


where 


An(a>) — 


ii - /— Y 

P n) 




and is the square 


root of the mechanical admittance for the nth bending mode 


0 n = tan 


jag, 

1 ■ (S) ! 


and 


'L, n'' 


(t) = / c L (|,t) h n (|) d£ 

v J * 


The total deflection at the vehicle nose is 


10 . 


y(0,t) = qRL 


An(“) e n 


n=l 


M rfV 


The power spectrum of the nose amplitude is 

<D y (0,a>) = lim 7^ F<{y(0,t)j> F*^y(0,t)} 


T— > CO 


2itT 


(*> 


(5) 
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where F^y(0,tj} is the Fourier transform of the nose deflection and the asterisk 
indicates the complex conjugate of the Fourier transform. By employing equa- 
tions (4) and ( 5 ), the spectrum of the nose displacement is 


00 00 

v<o>“>> - i 2 r 2 l 2 V y (ai) oo _ 

The function C m ^ n (a>) is a correlation function and is defined by 

Cm > n(a>) = 2^T *{*!., n^)} F *{%n^)) 


( 6 ) 


For a system with small damping and reasonably well separated natural frequencies^ 
the cross terms (m ^ n) in equation ( 6 ) may be disregarded since all contributions 
to the total response are small except in the neighborhood of the resonant fre- 
quencies. Thus, 


$y(0,O)) = q% 2 L 2 


> A n 

l— M n ^Dn^ 


c L,n(“>) 


(7) 


where C L ^ n (cn) = C n ^ n (cn). The mean-square nose deflection is 


— n 00 

y( 0 ,t) = J $y( 0 ,oj) cko 


(8) 


n^ 0 ^ reasonably constant in the vicinity of cDq and the damping is 

small, a satisfactory approximation to the mean- square nose deflection may be 
made. This approximation is 


y(o,t ) 2 


q 2 R 2 L 2 


00 

z 

n=l 


^n 




7n 


(9) 


In buffeting studies on elastic structures, usually the bending moment at 
some point on the structure, or the bending -moment distribution, is desired rather 
than the deflection of the structure. By using equation ( 9 ) and a set of coef- 
^^-Cients which relate the bending moment in the nth bending mode to the nose 
amplitude in that mode, an expression for the mean- square bending moment may be 


2k 


obtained. The bending moment at some 
deflection in the nth mode is 


longitudinal station £ 0 per unit nose 


^L 2 f l ° m(|) (| 0 - t) Ml) as 

J o v 

The expression for the total mean- square bending moment is 



i 2 


00 

2-n2 r 2 / 

L 2 / 0 m(l) - l) h n U) d| 

J n ' 

a l / 

1 

1 


n=l 



c L,n(“n) 


( 10 ) 


The term enclosed in the brackets in equation (10) 
and is abbreviated as Z n . The contribution of the 


is an effective moment arm 
nth mode to the total moment 


is 


a n 2 (^o) ~ l n 




k( r-) 

b cr/n 


q 2 R 2 L 2 C^ n (^) 


( 11 ) 


The right-hand side of equation (ll) is the general term of the series in equa- 
tion (10). The contribution of a single mode is essentially f, S ^^ d S + Q U a 
that would be obtained for a single-degree-of-freedom system subjected 
random loading. (See ref. 1.) Consequently, the total mean- square bending 
moment is a superposition of single-degree-of-freedom solutions, each mode being 
independently treated as a separate system. The right-hand side ° f 2 eq ^ 
may be conveniently separated into three parts. The first tern l n , as previ 

ously mentioned, is the square of an effective moment arm. ^ se ^ te ^ 
enclosed in brackets, is an admittance -type term. In particular, it is / 
times the maximum value of the mechanical admittance in the nth mode multiplied 


U line Id UXIC iucLA.XiULiiu — / C \ 

by its width at the one-half power point. The damping ratio (c cr ) n 


which 


appears in “this term is 
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where 


£a_ 

C cr 

£c_ 

Ccr 

£s_ 

Ccr 


aerodynamic damping ratio 


control-system damping ratio 


structural damping ratio 


In the notation of reference 7 the aerodynamic damping ratio is related to the 

effective aerodynamic damping derivative O by 

h J 

— - — ( 12 ) 
c cr 2 pk 

The remaining terms in equation (ll) are associated with the random aerodynamic 
loading. The function C L ^ n (cu) is a correlation function. Or more precisely, 

it is the correlation function of the random section lift coefficients for the 
nth mode. Although not mathematically exact, a convenient way of thinking of 
this function is that it is the power spectrum of an effective random aerodynamic 
coefficient. The correlation function should be more properly expressed as a 
function of reduced frequency rather than as a function of frequency. From 
dimensional considerations. 


where 


c L,n(“) - | C L , n ( k n) 




V 


(13) 


On substituting expressions (12) and (13) into equation (ll) the mean- square 
bending moment in the nth mode becomes 
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The variation of the mean- square 
a, given Mach number is of interest • 


bending moment with dynamic pressure for 
For the limiting case 



= 0 



= 0 



« q. 


anH for vanishing aerodynamic damping 




0 



In cases where both aerodynamic and structural damping must be considered, no 
simple proportionality relationship exists. 


Equation (15) is useful in scaling bending-moment measurements on dynamically 
scaled aeroelastic models to full-scale values. Since k n and L,n\ n ) 
be the same for both the model and the full-scale vehicle, the full-scale moment 
for the nth mode is related to the corresponding model value by 



where the subscripts M and F refer to model and full-scale configurations, 
respectively, and 

a F _ L F 2R F 2 ^n,F 2 
% LM 2 R M 2 Z n ,M 2 

For a dynamically scaled aeroelastic model, the ratio A F /A F is merely the 
reciprocal of the geometric scale factor Lm/W raised to the sixth power. 
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The total full-scale mean-square bending moment in terms of model bending 
moment is ’ 



It should be pointed out that expressions similar to equation (17) could be 
obtained for any quantity which is proportional to the displacement . 
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